VAND-TH-93-4 



2 



HIGGS-BOSON DECAY TO FOUR FERMIONS 

INCLUDING A SINGLE TOP 
QUARK BELOW THRESHOLD 

Heath Pois^'^, Thomas J. Weiler^, and Tzu Chiang Yuan^* 

^ Center For Theoretical Physics, Department of Physics, 
Texas A&M University, College Station, TX 77843-4242 

Astroparticle Physics Group, Houston Advanced Research Center (HARC) 
The Woodlands, TX 77381, USA 

^ Department of Physics and Astronomy, 
Vanderbilt University, Nashville, TN 37235 

^ Department of Physics and Astronomy, 
Northwestern University, Evanston, IL 60208 



Email address: pois@tamphys, weilertj@vuctrvax, nuhep::yuant 

* Present address: Department of Physics, University of California at Davis, CA 95616 



ABSTRACT 



The rare decay modes Higgs — > four light fermions, and Higgs — > single top- 
quark + three light fermions for mt < Mjj < 'imt, are presented, and phenomeno- 
logically interpreted. The angular correlation between fermion planes is presented 
as a test of the spin and intrinsic parity of the Higgs particle. In Higgs decay to 
single top, two tree-level graphs contribute in the standard model (SM); one cou- 
ples the Higgs to W+W^{r^ gMw), and one to gtopyukawa = mt/246GeV). 
The large Yukawa coupling for mt > 100 GeV makes the second amplitude com- 
petitive or dominant for most Mj/, mt values. Thus the Higgs decay rate to 
single top directly probes the SM universal mechanism generating both gauge 
boson and fermion masses, and offers a means to infer the Higgs-tf Yukawa 
coupling when — > is kinematically disallowed. We find that the modes 
pp Xtt{H thW^*^) at the SSC, and e+e" ^ ZorvV+{H ^ thW^*^) at fu- 
ture high energy, high luminosity colliders, may be measureable if 2mt is not too 
far above Mu- We classify non-standard Higgscs as gaugeo-phobic, fermio-phobic 
or fermio-philic, and discuss the Higgs— > single top rates for these classes. 



1. INTRODUCTION 



The existence of a Higgs scalar particle, or the onset of new physics at a 
few TeV or less, is a virtually guaranteed happenstance in particle physics [1]. 
Either possibility provides a foundation for the electroweak unification and the 
generation of masses. Of the two, the simpler possibility is the existence of 
the Higgs particle. Accordingly, a tremendous amount of energy has been, and 
continues to be, devoted to theoretical and experimental searches for signatures 
of the Higgs boson [2] . The four experiments at LEP have recently placed a lower 
bound of Mh ~ 57GeV@95%CL for a Standard Model (SM) Higgs [3]. With 
regard to the top quark mass, a direct SM lower bound of mt ~ 89 GeV has also 
been obtained from the CDF experiment [4]. 

Indirect SM upper bounds for mt and Mjj can be predicted by the theory 
based on quantum loop phenomenology. It is well known that heavy top quark 
loop corrections to certain low energy and electroweak (EW) observables (for 
example, the p parameter) are proportional to ml and thus the quantum effects 
are quite sensitive to mt. The SM consistency of all the low energy experimental 
data requires mt < 182 GeV@95%CL, with a central value of mt = 125 ±30 GeV 
[5]. On the other hand, the mass dependence of a heavy Higgs loop correction 
varies as InM^ in the SM. This is the famous one-loop "screening rule" first 
recognized by Veltman [6] . Since the dependence of quantum loop effects on the 
heavy Higgs is only logarithmic, low energy observables are not very sensitive to 
Mjj. However, recent analysis [7] indicates that a weak upper bound {Mh ~ 
300 — 500 GeV) for the SM Higgs can be deduced from current low energy 
experimental data which is suggestive of a 'light' Higgs. In models with a broken 
supersymmetry, the tree level mass of the lightest Higgs generally lies below the 
Z mass [2]. However, the large radiative corrections to the Higgs boson masses 
due to a heavy top quark [8] may alter this situation significantly. 

Finding a signature for an 'intermediate mass' Higgs (Mw < Mjj < 2Mz) 
is particularly problematical. It appears that rare decay modes of the Higgs 
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such as — > Z7, 77, t''"t~, 66, ©7 offer the greatest promise for providing an 
experimental signature of the particle's existence [9] . The decay of a Higgs with 
mass above the two W threshold is dominantly to two VF's, and so one might 
believe that a Higgs with mass not too far below twice the W mass may have 
a significant branching fraction to decay through WW* or intermediate 
states {W* denotes a virtual W). 

Marciano and Keung [10] have shown that for a standard model Higgs with 
mass below but near the two W threshold, the branching ratio for the decay 
H — > VFVF* becomes significant; in fact, this decay mode dominates all others if 
the decay to tt is kinematically forbidden (i.e. Mb: < 2mt). Barger et. a/. [11] 
have considered the — > four fermions decay mode in the massless fermion 
limit, and have made a thorough study of the SSC signal and backgrounds. Their 
conclusion is that only if mt, Mjj > 150 GeV will the decay signal be seen 

above the tt — > WbWb and W continuum background channels. 

Our results here are complementary to the work in [11], and extend that work 
by including a single massive top quark in the final state. As in [11], we allow 
both W^s to be virtual, and let phase space optimize the sharing of virtuality 
between the two W^s. Our exact results for Higgs four fermions allow us to 
exhibit the correlations among final fermion energies and angles that test the 
presumably scalar nature of the parent Higgs particle. As an example of final 
state correlations, we display the dependence of the decay on the angle between 
the decay planes defined for each vector boson. The analogous azimuthal angular 
dependence for the process pion e+e^e+e^ was calculated over thirty years 
ago [12], and ultimately provided the signature distinguishing between a scalar 
and a pseudoscalar 'pion' [13]. 

In Sec. 2 we present our result for the massless fermion limit (which agrees 
with the result in [11]), and we discuss a class of nonstandard Higgses which do 
not couple directly to fermions. Such 'fermiophobic' Higgs decay through loops to 
two fermions, or through two gauge bosons to four fermions. The 'fermiophobic' 
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Higgs branching ratio to four fermions is large, even if mass kinematics require 
both gauge bosons to be virtual; a calculation of this branching ratio requires 
the exact formulae presented in this paper. In Sec. 3 we define the asymmetry- 
parameters relevant to the angular correlations between the fermion planes, and 
explore the Mh and virtual gauge-boson mass dependence of these parameters. 

In Sec. 4 we exhibit the rate for the tree level decay of the Higgs to a single 
heavy quark (e.g. top) plus three light fermions. Since the Higgs particle has a 
direct Yukawa coupling to fermions which scales with the fermion mass, this latter 
calculation includes a second graph with a potentially large contribution to the 
amplitude, and therefore provides an important contrast to the massless fermion 
case. We find that due to the large Yukawa coupling for large m^, the branching 
ratio to single top (e.g. H tbsc) can be almost competitive with the dominant 
massless fermion decay mode H — > WW 4/, and is certainly competitive 
with the rare H — > 77, Z7, bb, t~^t~ decay modes for a range of Mjj, mt values. 
In Sec. 5, we present a brief discussion of the signature and backgrounds for 
the H tbW* mode. We argue that the signal may be detectible at hadron 
colliders for Higgses produced in association with a tt pair, and at future high 
energy, high luminosity e'^e~ colliders. Interesting non-standard Higgs models 
are classified in Sec. 6, and the branching ratios to single top of these baroque 
Higgses are discussed. We conclude in Sec. 7. Detailed formulae for the decay of 
a Higgs boson to four fermions, including the massive heavy top, are collected in 
an appendix. 



2. HIGGS DECAY TO FOUR MASSLESS FERMIONS 

For H WW, ZZ 4/ (where the W^s and Z's may be real or virtual), the 
relevant Feynman diagram when mj = is shown in Fig. 1(a); Fig. 1(b) obviously 
does not contribute in the massless limit to the SM Higgs due to the vanishing 
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Yukawa coupling. (See however Sec. 6 for non-standard Higgses). The result for 
the Higgs width to four fermions is presented in the appendix. The differential 
width is symmetric in the invariant masses and ^^e W or Z gauge 

bosons, and we show the gauge boson mass spectrum in Fig. 2 for the process 
H WW udsc + udsc. One can observe that once Mjj > M^, a bimodal 
distribution results. The twin peaks correspond to one or the other W being 
on-shell; for invariant mass values Minv between the peaks, kinematics requires 
both Ws to be off-shell. It is clear that the relative contribution with both W^s 
off-shell is small for an intermediate mass Higgs. This observation is borne out 
in a comparison between the exact four fermion width, and the one gauge boson 
on-shell (OGBOS) approximate width [10], which is obtained by replacing one of 
the Breit-Wigner factors ((Q^ - M^f + M^F^)"! in the decay width with the 
narrow width approximate (NWA) form 7r5{Q^ — My) / MyTy ■ In agreement with 
[11], we find that putting one W on-shell is a good approximation to the total 
decay rate for Mjj > 100 GeV. For Higgs decay through two Z's, we find that 
putting one Z on-shell approximates well the true decay rate if Mfj > 115 GeV. 
In conclusion, we find that the rate for a SM Higgs particle to decay to four 
fermions is a percent of that for decay to bb, for Mjj ~ 100 GeV, and falls 
rapidly for lighter Higgs masses; for Mh > 100 GeV, the OGBOS approximation 
gives an accurate rate. 

A bit of caution is required when applying the NWA to the Breit-Wigner inte- 
grals of Eqn.(A.21). If one Breit-Wigner is replaced by the NWA delta function, 
a multiplicative factor of two must be introduced to account for the probability 
for either W to be on-shell. Above the two W threshold, this factor of two must 
be removed since one W on-shell no longer precludes the possibility of the second 
W also being on-shell. If Eqn.(A.21) is used without approximation, there is no 
need to concern oneself with this extra counting. 

The four-fermion mode becomes particularly interesting in models where the 
two-fermion modes are suppressed at tree level. Recall that in the SM with three 
families of quarks and leptons, the Higgs doublet plays a double role of gener- 
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ating masses for both the gauge sector and fermion sector. In fact, there is a 
third role: the vacuum expectation value (VEV) of the Higgs field breaks not 
only the SU{2) x U{1) gauge symmetry but also the flavor symmetry [t/(3)]^ 
(the three copies of the five fields II, Ql, ^r, ur, and da). The Higgs mechanism, 
originally employed to give nonzero masses to the gauge bosons while maintain- 
ing a renormalizable theory, may not be the source of the fermion masses or the 
flavor symmetry breaking. The independent generation of fermion masses and 
gauge bosons masses is conceptually possible [14]. A suppression of the Higgs 
couplings to fermions is easily enforced by introducing a discrete (f) — > —(f) sym- 
metry; the symmetry forbids Yukawa couplings [15]. We call the resulting Higgs 
fields 'fermiophobic', for although the fields couple to the gauge bosons at tree 
level, their coupling to fermions occurs only through loops and scalar field mixing 
(The loop and mixing-induced couplings are allowed since < (f) >^ breaks the 
discrete symmetry). 

It is a logical possibility then that the sole Higgs particle introduced in the 
standard model is itself 'fermiophobic'. If so, the rate for a Higgs of intermediate 
mass to decay to four fermions via two virtual gauge bosons (Fig. 1(a)) competes 
favorably with all other decays, namely Higgs two fermions through scalar 
mixing and/or loop graphs, Higgs — > 77, or 7Z through a, W or charged Higgs 
loop, and Higgs — > gg through one-loop and mixing, or through two loops. In 
fact, the decay rate to four fermions exceeds the rate to 77 or 7Z if the Higgs 
mass exceeds 75 GeV. Branching ratios for a fermiophobic Higgs may be inferred 
from the SM branching ratios presented in Fig. 2 of ref. [11]. If an intermediate 
mass Higgs is discovered, a comparison of its two fermion decay rate to its four 
fermion or two photon rate will immediately provide important information on 
the mechanism(s) of mass generation. [16] 
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3. ANGULAR CORRELATIONS BETWEEN THE FERMION PLANES 

Next we turn to the correlations among the final state of four massless 
fermions implied by the matrix element of Eqn.(A.4). The reader may recall 
that parity conservation implies that the two photons resulting from a scalar 
particle decay have parallel polarization directions, while the photons from a 
decaying pseudoscalar particle have perpendicular polarizations [17]. A similar 
result holds for decay of scalar and pseudoscalar bosons to two Z's or two 1^'s 
with momenta Qi and Q2 and polarization vectors e'^{X) and £2 (A'). The matrix 
element is linear in the two polarization vectors. 

In the case of even-parity Higgs field surviving spontaneous symmetry break- 
ing, the decay proceeds at the tree level with the two polarization vectors con- 
tracted ei(A) • €2 (A'), which tells us that the two gauge bosons emerge with the 
same linear polarization. In a helicity basis, the two gauge bosons emerge with 
the same helicity. For massive gauge bosons, there are of course three physical 
polarizations, and a fourth 'scalar' polarization for off-shell vector bosons which 
vanishes when coupled to a conserved current. In the case of the odd-parity 
neutral pseudoscalar decay, CP invariance of the gauge boson and scalar boson 
sectors is sufficient to ensure no tree level coupling [2] ; the process is loop induced 
with (pFF being the lowest dimensional operator, just as in the tt'^ 'j'j case. 
The Lorentz invariant resulting from this operator, linear in the two polarization 
vectors, is <^^waB'^i{^)'^^{X')QiQ2 (times a coefficient down by a compared to the 
tree level scalar case). In the pseudoscalar particle rest frame this invariant is 
proportional to ei(A) x €2 (A') • {Qi — Q2)i which tells us that only transverse polar- 
izations, or equivalently, left and right helicities, are produced. Furthermore, the 
linear polarization vectors are perpendicular to each other. In a helicity basis, 
the two helicities are again identical. We expect the final decay planes to remem- 
ber the vector boson polarizations, just as in electromagnetic decay of scalar or 
pseudoscalar particles to two e+e" pairs [12]. Since a scalar parent particle gives 
a Lorentz invariant combination of polarization vectors differing from that of a 
pseudoscalar parent, we expect the angular orientation of the ffnal fermion decay 
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planes to be an indicator of the parity of the parent particle. 

So far, our discussion of polarization has concerned itself only with vector 
polarization. In fact, since the polarization of the gauge bosons cannot be directly 
measured, a true description of the gauge boson polarization requires a density 
matrix; the vector polarization is just the diagonal part of the general density 
matrix. The density matrix for an off-shell vector boson is a 4 x 4 matrix. Since 
our calculation couples the gauge boson to a conserved final state current (in 
the massless fermion approximation), only the 'physical' 3x3 sector contributes. 
When either gauge boson is put on-shell and treated as an observable particle, as 
with the OGBOS approximation, the correlations contained in the off-diagonal 
elements are lost! Of course, all of the correlations implied by the density matrix 
formalism are contained in the exact result of Eqns.(A.17-21), which we use. 

Let us define the azimuthal asymmetry parameters ai and 0:2 according to 
the following formulae: 



dr _ T 

'd4>~ 2tt 



1 + aiiy) COS0 + a2iy) cos 20 



(1) 



(f) is the angle between the two planes, where each plane is defined by the two 
fermion momentum vectors resulting from the decay of a virtual or real gauge 
boson, as measured in the four fermion CMS (i.e. the Higgs boson rest frame), (f) 
is invariant under boosts in the direction of the momentum of either gauge boson. 
V stands for either the WW decay mode, or the ZZ mode. The expressions for 
a\{V) and a2{V) for scalar Higgs decay are given in Eqns.(A.17,A.18). a\{V) is a 
parity violating asymmetry and a2{y) is parity conserving. We show in Fig.3 the 
asymmetry parameters as a function of the Higgs mass Mjj^ for the decays H — > 
ZZ — > ddss, and H — > WW — > udsc+udsc. For convenience, we have chosen final 
states such that identical fermions, requiring antisymmetrization, do not occur, 
and such that interference between the two-W graph and the two-Z graph does 
not occur. It is seen that the asymmetries are largest for the intermediate mass 
Higgs, and peak near the threshold values Mh — 2Mw, '2Mz for q;i,2(VF), ai^2{Z) 
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respectively. At threshold the two F's are nearly on-shell and at rest, so the 
three nonzero helicity amplitudes are equal and the off-diagonal density matrix 
elements are maximized. In particular, the parity violating asymmetry a-\_{V) can 
be of for the intermediate mass Higgs, and may provide a clean consistency 
check on the scalar nature of the decaying Higgs particle. In comparison, from 
Fig. 3 we see that the parity conserving asymmetry is reduced by an order of 
magnitude or more. 

Unfortunately, the large asymmetry cti is difficult to determine experimen- 
tally. The azimuthal angle may be written as cos0 = ni -712, where ni^2 Sire the 
normals to the two WW or ZZ decay planes. To determine these normals, an 
experimenter must have charge and possibly flavour identification of the emitted 
leptons and/or jets. An f ^ f ambiguity results in a n <-> —n or (p ^ tt — (p am- 
biguity, and washes out the cos0 asymmetry ai. On the other hand, the cos2(/) 
term is even under this ambiguity, but a2 is small except near VV threshold. 

The asymmetries have been calculated before for on-shell W^s and Z's [18], 
and for the QCD background process gg, qq — > ZZ — > 4/ [19]. Applying the NWA 
to our more general results, Eqns. (A. 17- A. 18), or referring to the literature, one 
finds simple expressions: 



CKi on — shell V s) = — — - — ^ n.^ — a i 

^ ' 16 \Cr + Cl) [(M^ - 2M2 )2 + 8M4] 



(2) 



2M^ 

At threshold a2{y) = 1/6 (evident in our figure); above threshold a2{W) and 
a2{Z) fall rapidly (asymptotically, like {Mv/Mjj)^) to less than two percent at 
Mff > 300 GeV. What is new in our calculations are the below threshold results: 
a2{V) falls as Mjj moves below the VV threshold, but slowly. When kinematics 
require one (two) V off-shell, the asymmetry is ~ 10% (7%). For ai{V), the 
threshold value is ^ ^ cl+c'l ) • are given in the appendix. For the 
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V = W case, (where parity violation is maximal, and Cl = 1, Cr = 0) ai{W) = 
= —0.925, nearly providing a zero in ^ at cos0 = 1. Extremization of 
aiiy) with respect to Mh reveals that |Q!i(y)| achieves its maximum value a bit 
above threshold, at Mh = 2.35My, i.e. at 188(214) GeV for the W{Z) boson. 
This slight displacement is also evident in our figure. The above threshold fall- 
off is not so rapid as for going asymptotically like {My/Mjj)'^- The below 
threshold fall-off for cti is also weaker than for 0:2, dropping by 15%(30%) when 
one (two) Vs are kinematically moved off-shell. 

The asymmetries arise from interference among the possible H 
helicity amplitudes. In terms of helicity amplitudes, the dependence of ai and 
a2 on Mjj is qualitatively easy to understand. Only interferences between ±1 
helicities and zero helicities contribute to cti, while only interferences between 
±1 helicities and ±1 helicities contribute to Q!2- At threshold, all helicity states 
are equally populated and the relative interferences are maximized. For an off- 
shell V, the longitudinal mode is slightly more populated than the transverse 
mode, reducing the relative interference and thereby, the asymmetries. Above 
threshold, the transverse modes are greatly suppressed relative to the longitudinal 
mode, asymptotically like {My /MhY , and so the asymmetries ai and 012 fall like 
{Mv/Mnf and (Mv/Mh)^, respectively. 

As may be seen from Eqns.(Al7-20), the asymmetry coefficients are indepen- 
dent of the final state except for ai{Z). ai(Z) is maximized when Z*Z* two 
down-type quark pairs. For each up-type quark pair in the final state ai{Z) is 
reduced by 0.73; for each charged lepton pair by a factor of 0.17. In fact it is easy 
to see that ai{V) is maximized when the Vff coupling is purely chiral (either 
left-handed or right-handed) so that parity violation is maximal, and goes to zero 
for a coupling that is pure vector or axial vector. 

The nearness of sin^ 9^ to the "magic" value 1 /4 means that the vector cou- 
pling of the Z to a charged lepton pair is nearly zero; this in turn means that 
the asymmetry ai{Z) with a in the final state is nearly zero. This circum- 
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stance is most unkind, since the Z — > mode has the best signature, both 
in terms of beating backgrounds, and beating the ±n ambiguity. Prospects for 
reconstruction of the Higgs and its decay planes are better at e+e" machines 
than at hadron coUiders. 

As explained earlier, a CP-odd pseudoscalar meson does not couple to the 
longitudinal V%. Accordingly, it has a vanishing cti, but a maximal and negative 
a2 value of ^ [18]. Such would be the asymmetries for a technipion, for example. 

Incidentally, it is just this angular orientation of decay planes that was used to 
definitely determine the parity of the neutral pion thirty years ago [13]. Perhaps 
after extraordinary effort physics history will repeat itself. Or perhaps other 
observables, such as the fermion-antifermion lab energy asymmetries Ej^/E^ or 
averaged rest frame polar angles [20], will prove more useful. 



4. HIGGS DECAY TO FOUR FERMIONS INCLUDING SINGLE TOP 

The coupling of the standard model Higgs to a fermion scales as the fermion 
mass over the electroweak VEV, ~ 250 GeV. Therefore, unless a fermion is 
heavy on the scale of the W or Z, its coupling to the Higgs is negligible. Here 
we consider the modification of our decay rate to four fermions when fermion 
masses are included. We take one fermion to be massive, and continue to neglect 
the other fermion masses [21]. The relevant example we consider is the decay of 
a Higgs with mass above single top threshold but below ti threshold. The decay 
mode is — > top (mass mt) plus three light (on the scale of the W boson, the 
6— quark qualifies) fermions. In models where the Higgs is generated dynamically 
as a scalar tt bound state [22], or where the top and Higgs masses are tuned to 
cancel divergences [23], the Higgs mass is expected to lie in this range [m^, 2mt\. 
We also note that the branching ratio of toponium to a Higgs in this mass range 
is expected to be a few percent via the channel. However, toponium is not 
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expected to exist if mt exceeds about 140 Gev, since then the top quark hfetime 
is shorter than the bound state formation time. Thus, if the tt threshold exceeds 
the Higgs mass, then H single top may offer the best hope for extracting the 
Htt coupling. 

In addition to the graph of Fig. 1(a) already considered, there is the additional 
graph of Fig. 1(b), proportional to the Htt coupling. Since Z bosons conserve fla- 
vor and therefore cannot produce a single top, the only gauge boson contributing 
in Fig. 1(a) is the W. For improved accuracy, we keep the mass of the 6-quark in 
the phase space, although it is omitted in the matrix element. 

The expression for the width is presented in Eqns.(A.l)-(A.6). The width for 
the process H WW tbsc + tbsc for various values of mt/Mi{ is displayed 
in Fig. 4. Single but not double production of top requires 0.5 < mt/Mu < 1.0, 
and various choices are indicated in the figure. From direct CDF top mass limits 
[4] and from radiative correction limits [5], 89GeV < mt < 182 GeV at 95%CL. 
In our analysis, we relax the CL and assume that mt < 250 GeV absolutely. For 
comparison, we include the massless limit {mt/Mff = 0) and the H — > bb) width 
in Fig. 4. Were the Higgs mass above tt threshold, the H ^ tt width would 
be the bb curve scaled up by /3'^/^(mt/mfe)^, where (3 = ^1 — 4mf/Mjj. The 
massless limit describes, for example, the process H — > WW udsc + udsc, 
and is roughly 25% of the H WW width. One can see that for mt/Mfj ~ 0.5 
(precisely at ti threshold) the partial width for H tbsc + tbsc is reduced by 
only a factor of 10 or less compared to the massless limit for Mjj ~ 200 — 500 
GeV. For larger values of mt/Mjj, the partial width falls off rapidly. The ratio 
n^t/Mff determines how far off-shell the virtual top is, and so determines the 
qualitative features of the rate. Fig.5 amply demonstrates this strong sensitivity. 
Starting at mt/Mff = 0.4 {Mh = 400 GeV), the rate for H tbsc + tbsc drops 
by C( 10^) at mt/Mu = 0.55. For decreasing values of Mjj and fixed mt/M^, the 
drop is increasingly more severe due to diminishing phase space. If mt/M^ = 0.6 
and Mh > 200 GeV, T{H tbsc + tbsc) /r{H ^ r+r") = C»(l), which implies 
a drastic reduction, since H — > t~^t~ is well known to be a rare decay mode. 
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An important question to consider at this point is the relative magnitude 
of the couphngs that contribute to the partial width. For large mt, one would 
guess that the Yukawa contribution should be significant. This relationship is 
indeed demonstrated in Fig. 6; the three curves show the separate partial width 
contributions arising from the Yukawa (Y), gauge (G), and Yukawa + gauge (1^ + 
G) terms including the interference. The two graphs of Fig. 1 are individually 
gauge invariant so the separation into Y and G is meaningful. In Fig. 6, mt/M}j is 
chosen to be 0.6, so 89 GeV <mt< 250 GeV confines 150 GeV < Mh < 420 GeV 
as the physically relevant region. One can observe that Fy is the dominant 
contribution for this range of Mjj by an order of magnitude over Vq. At fixed 
Ty will dominate even more for larger values of mt/Mu- 

The structure of the curves can be understood as follows: In the gauge 
amplitude (Fig. la) the intermediate W connecting to {tb) is kinematically al- 
lowed to venture on-shell when Mh > Mw > mt + mi,. In Fig. 6 this occurs for 
Mw < Mff < 125 GeV, and we see structure in Tq in this region. In the gauge 
and Yukawa amplitudes the intermediate W connecting the massless qq pair is 
kinematically allowed on-shell when Mh > Mw + mt + mf). In Fig. 6 this occurs 
for Mjj > 215 GeV, and we see a sharp rise in both amplitudes at this value. 



5. SIGNATURE AND BACKGROUND 

The final states for H tbW-^*^ bW+^%W-^*^ will be the same as those 
from the QCD process gg ^ tt ^ bW+^%W-^*\ thus making the H t3f 
mode difficult to extract at a hadron collider. The signal to background issue 
here is quite similar to the well known difficulties in extracting an above-threshold 
H ^ tt signal from a QCD background. In fact, within the top mass resolution 
expected at LHC and SSC detectors [24], the bW~^*^ accompanying the single top 
may well mimic an associated i. With a bW^*'> / associated-top misidentification. 
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the signature difficulties are identical for the H — > single top and the H ^ tt 
possibilities. However, unlike the H ^ tt branching ratio, the H single top 
branching ratio is small, making signal extraction much more difficult. 

What may be fruitful at a hadron collider is to trigger on ti produced in 
association with Higgses. The dominant Higgs inclusive production mechanism 
below Mh = 500 GeV is gluon-gluon fusion, which yields a cross section of ~ 40 
pb at SSC energies {^/s = 40 TeV), fairly fiat for Mh < 500 GeV. The Hti cross 
section at SSC energies, also dominated by the gluon-gluon fusion mechanism, 
is ~ 4 pb at Mh = 150 GeV, falling to ~ 0.3 pb at Mh = 500 GeV; it is rather 
independent of the top mass in the range of our interests. Thus the loss in event 
rate is a factor of only ~ 10 to 130; on the order of 10^ Htt events per year may 
be expected assuming the standard SSC integrated luminosity of 10^ pb^V^SC 
yr (10^ sec). The Htt cross section at the LHC {^/s = 16 TeV) is down from that 
of the SSC by an order of magnitude or more (a function of masses). However, 
the higher design luminosity at the LHC may argue for a number of Htt events 
comparable to the SSC. It has been shown recently [25] that triggering on an 
isolated secondary lepton, e or from decay of the associated t or t, renders 
discovery of the intermediate mass Higgs feasible at the SSC, via the rare decay 
mode H — > 77, {BR ~ 10~^); in an analogous fashion, one may hope that 
extraction of the rare decay mode presented here, H — > tbW^*\ may also be 
feasible in association with a.tt ^ isolated trigger; a, pp ^ tt{H tbW^*^)X 
event (or above tt-threshold, the pp tt{H tt)X event) would lead to a 
spectacular number of high pT jets and Icptons. If the Higgs is not reconstruct ible 
in the single top mode, then an excess of spectacular events may be the only 
signal from the mode. Possibly helpful for signal enhancement is the fact that 
W's from the Higgs decay chain will be predominantly longitudinally polarized 
[26], whereas "background" W's produced by light quarks will be transversely 
polarized. 

The order pp titi QCD background to our process has been calculated 
recently [27]. The result is crgcDipp ^ ttttX) = 1, 0.4, 0.2, and 0.07 pb at 
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the SSC for mt = 110, 140, 170, 200 GeV; and 0.1, 0.03, 0.008, and 0.003 
pb for the same mt values, at the LHC. Thus, at either the SSC or the LHC, 
a{pp — > ttH X) / aqcoipp — ^ ttttX) ~ 0.3 to 50, depending on the Higgs and top 
masses, in the range of our interest (100 GeV < Mh < 500 GeV, 90 GeV < 
mt < 250 GeV). Thus a BR> 10"^ in the H tbW^*^ mode may in principle 
be observable. From figure 4, if M}j > 160 GeV, then mt < 0.55M}j appears 
to be needed to obtain such a branching ratio. A concern is that there is no 
obvious signature distinguishing the from H tbW^*\ from the second tt 
in the QCD signal; all tops will tend to be produced near threshold. The extra 
combinatoric possibilities available with four tops is also problematic. A Monte 
Carlo simulation is required to quantitatively determine detection capabilities for 
the H tbW'^*^ mode. 

At a future e^e~ collider, it may prove possible to probe a smaller BR{H 
tbW), and hence a larger mt/Mu ratio, if the integrated luminosity is sufficiently 
high, and if the ti and titi backgrounds are more controllable. The associated 
Z+H and/or WW fusion cross section for production of a 100 GeV SM Higgs is 
0.30, 0.14, and 0.24 pb for the values = 200, 500, and 1 TeV, while for a 200 
GeV Higgs it is roughly 0.1 pb over the range ^/s = 500 GeV to ^/s = 1 TeV, 
and for a 500 GeV Higgs it is 0.03 pb at = 1 TeV [2]. Thus, with 10/6" 
of luminosity, 300 to 3000 Higgs events per year are anticipated. A 100/6~^/yr 
machine has even been discussed, for which the rates are larger by a factor of 
ten. We conclude that branching ratios as small as 10^"^ may be detectible. 

Finally, we comment that if Mjj exceeds tt threshold, as opposed to the mass 
range mt < Mjj < 2mt considered in this paper, then 

(7qcd{pP ^ ttttX) ^ ) \ )^ \ ) 

which suggests that the Hti coupling may be experimentally accessible if BR{H 
tt) exceeds 10"^. Using SM formulae, one finds BR{H ^ tt) ^ T{H 
tt)/T{H ti,ZZ,W+W-) = r(l -r)i/[l + r(l - r)i], where r = Amf/Mjj, 
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and terms of order Mz/Mh have been neglected. This branching ratio can be as 
large as 30% (at r = 0.4, or mt ~ M^/S), clearly warranting further study of this 
mode. As with the sub-threshold process H ^ tX, a Monte Carlo simulation is 
required to determine signal viability. 



6. FERMIOPHOBIA, GAUGEOPHOBIA, AND FERMIOPHILIA 

Given the possibility that gauge boson masses and fermion masses may 
arise from different mechanisms, one can imagine scalar fields (i) with cou- 
plings to fermions suppressed relative to the SM values ruf/vEWi (h) with cou- 
plings to WW and ZZ suppressed relative to the SM values Q^vew/'^i ^■nd 
(p'VEW 1^ cos^ 6*^), respectively, and (iii) with enhanced couplings to the fermions. 
We will label these scalars fermiophobic, gaugeophobic, and fermiophilic. One 
can easily invent models of each kind, and find examples of each kind in the 
literature. We have already mentioned the (/> — > — fermiophobic scalar. Any 
scalar in an SU{2) representation other than doublet will also be fermiophobic, 
since the SU{2) invariance then forbids the (j)^ r operator. 

The (t)^WW,(l)^ZZ,{(j)*^(j))W, and {(t)*^ (fyZ couplings are fixed by the 
gauge coupling g. However, the (f)WW and (pZZ couplings are proportional 
to < (f) >. Accordingly, the latter are suppressed if \/2 < > is less than 
VEW = 246 GeV. In fact, the measured value of the p-parameter tells us that 
VEVs of representations larger than doublet are much smaller than veWi so 
scalars from representations larger than doublet are gaugeophobic as well as 
fermiophobic. Higgs doublets may themselves be gaugeophobic, since in a multi- 
doublet model, v'^y^ = X^doublets ^l) sums the positive contribution from each 
doublet. Further models abound. As mentioned earlier, in models which con- 
serve CP in the Higgs and gauge-boson sectors, the pseudoscalar Higgses are 
gaugeophobic [2]. In general two-doublet models, the charged Higgs and the 
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single CP-odd Higgs are gaugeophobic [2] . 

Fermiophilic scalars may result from multi-scalar mixing. A well known 
example is the general two-doublet model, in which some Yukawa couplings 
are enhanced by VEV ratios, while others are suppressed (fermiophobia again) 
by the inverse ratio. (Particularly well motivated two-doublet models are the 
minimal supersymmetric model invented to stabilize the weak scale, and the 
Peccei-Quinn model [28], invented to solve the strong CP problem.) There is 
also the possibility of scalar fields in representations which do not get a VEV. 
These non-Higgs scalars would have no tree-level (f)WW or (f)ZZ couplings (i.e. 
gaugeophobic) and arbitrary Yukawa couplings (fermiophobic or fermiophilic). 
Such "extra" scalars have been invoked in many theoretical contexts [29]. 

Finally, we mention hybrid models where the VEV of a fundamental scalar 
field generates one mass scale, while other masses arise from a condensate or 
from mixing with heavy singlet fermions [30]. Such models incorporate the best 
attributes of technicolor or fermion singlet masses with the simple attributes of 
Higgs physics, but with an increase in field content and a decrease in aesthetics. 
What is clear is that when a Higgs particle is finally discovered, its branching 
ratios will be very revealing. From the SM branching ratios of Fig. 2 in reference 
[11], one may infer the dominant decay modes of non-standard Higgses: for 
a neutral fermiophobic Higgs, omit the two-fermion and two-gluon modes; for 
a neutral gaugeophobic Higgs, omit the four-fermion, WW, and ZZ modes. 
Ignoring the loop-and mixing-induced H ^ ff rate and the two-loop H ^ gg 
rate, one learns that the dominant decay mode for a fermiophobic Higgs is 77 
for Mh < 80 GeV, and four fermions for Mr > 80 GeV (via W*W* for 80 
GeV < Mh < 100 GeV, via W*W for 100 GeV < Mh < 2Mw, and via WW 
above the production threshold). Estimates of the rates for the induced // and 
gg modes require calculations within a specific model. For a heavy gaugeophobic 
Higgs, the tt mode is dominant over the WW mode. 

At the SSC machine with energy ^/s ~ 40 TeV, the gg — > top-loop —>■ H 
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chain is expected to be the dominant production mechanism for a standard Higgs 
with Mjj < 5mt (as here). However, a fermiophobic Higgs may in fact have 
as its dominant production mode WW fusion, or W* — > WH, with a smaller 
production rate. At an e'^e~ machine, the standard Higgs is expected to be 
produced via Z* ^ Z + H for up to AOOGeV + 0.6Mjy , and by WW fusion at 
higher [2]. A gaugeophobic Higgs would have a suppressed production rate, 
either in association with bb or tt (if kinematicaly allowed), or via a top- loop. 

Let us now discuss changes in the rate for H ^ t+ 3 fermions when the Higgs 
is non-standard. Let v/^/2 denote the VEV of the non-standard scalar multiplet. 
Then the HWW coupling is g'^v/2 and the Htt coupling is an arbitrary constant 
qy- If we make the reasonable assumption that all Higgs fields couple to top with 
the same sign, then mt = Xlscalars 9yv/V^ provides the constraint gy < \/2mt/v; 
we write gy ~ x/S^wt/^S < ^ < 1. li v < veWi the Higgs is gaugeophobic; 
if ^/v < 1/veWj the Higgs is fermiophobic; and ii ^/v > 1/veWj the Higgs is 
fermiophilic. If H is in fact a scalar without a VEV, then setting v equal to zero 
yields the appropriate result: decoupling from VF+VF" and an arbitrary coupling 
to tt. The rate for non-standard H ^ t+ 3 fermions in terms of v and ^ can be 
obtained from Fig. 6: Fy scales as (^vew/^)'^: Tg scales as (v/vew)'^: and the 
interference term scales simply as ^. For a truly gaugeophobic Higgs, Fy is the 
total width, while for a truly fermiophobic Higgs, F^ is_ the total width. 



7. SUMMARY 

We have presented the tree-decay of the SM Higgs to four fermions, including 
the possibility of a single massive fermion in the final state in order to explore 
the potentially large Yukawa contribution to the decay rate. In the massless 
limit, we have demonstrated that if both gauge bosons are off-shell, the rate is 
considerably reduced. If one gauge boson can go on-shell, it will, and therefore 
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the OGBOS approximation is valid for an intermediate mass Higgs. We have 
also explored how a measurement of the angular correlation of the decay planes 
of the final state fermions may lead to a determination of the intrinsic parity of 
the Higgs, even if Mh < 2Mw- 

Next, we have shown that for a Higgs mass below ti threshold, the decay 
rate to single top is generally dominated by the Yukawa coupling, and thus tests 
the SM mechanism for fermion mass generation. For a Higgs just below the tt 
threshold, the single top rate is only O{10~^) down from the massless mode where 
the gauge coupling is the only contribution. As mt/Mff grows, the overall rate is 
dramatically reduced. It is evident (see e.g. Fig.4) that the H — > tbW^*^ process 
remains competitive in rate with other potentially detectable rare decay modes 
such as H ^ 77, Z'j, bb, t^t~ over much of the allowed range of M^, nit/MH- 
At a hadron collider, the signal to background is enhanced by triggering on pp 
tlHX signatures, making H — * single top possibly observable if 2mt^l.l Mh', 
a Monte Carlo simulation is needed to quantitatively compare the signal to the 
0{ag) titi QCD background. (Prospects appear brighter for measuring T{H — > 
tt) via pp — > ttHX, should Mh exceed 2mt.) It appears that H — > single top 
may have a better chance for detection at a future high energy, high luminosity 
e+e" collider, where branching ratios as small as 10~^, i.e. 2mt^l.2 Mh, may 
be measureable. 

Finally, we have shown that the rate for Higgs — > single top is a sensitive 
measure of any "non-standardness" in either the Higgs-gauge boson and/or the 
Higgs-fermion sectors of the true theory. We havediscussed the implications of 
"non-standardness" in these sectors for the general issue of fermion and gauge 
boson mass generation. 
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APPENDIX 

The Lorentz invariant matrix element for particle decay to a four body final 
state depends on five independent variables. We find it convenient to use as 
variables the invariant mass of each fcrmion pair connected by a common fermion 
line (there are two of these), the polar angle of the particle momentum with 
respect to the pair momentum direction, evaluated in the pair center of mass 
frame (there are two of these) , and the azimuthal angle describing the orientation 
of one plane defined by paired fermion momenta relative to the plane defined 
by the other paired momenta (there is one of these). For the case with all 
massless final state particles, where only graph 1(a) contributes, the first four 
of these variables are just the virtual W or Z boson invariant masses, and the 
fermion momentum direction relative to its parent W or Z momentum direction, 
evaluated in the W or Z boson rest frame. The Lorentz invariant phase space in 
terms of these variables is 

M% {MH-^/Qir 2n 1 1 

J dlips^"^^ = J dQl J dQ\ j d(j) J dcos 9\ J dcos 02 

{mt+rrib)^ 0-1 -1 

Xy^iM]j,QlQl)Xy\Qlmlml) 



(Al) 



where A(a, 6, c) — -\- 1? + — 2{ab + be + ca) is the usual triangle function. 
(An alternate choice of the five independent phase space variables is discussed in 
[31]). The width is 



Y{H ^ {th + ih) + //) = ^ / dlips^^^ 



spin 



{A.2) 



where 



(A.3) 



MaMl + MbMl + 2ReMaMl 

spm spin '- ^ 

expresses the summation over the squares of the amplitudes from Figs. 1(a) and 
1(b), and the interference of the two amplitudes. The expressions for these 
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squared matrix elements are 



J^MaMl = rg^M^Bw{Ql)Bw{Ql) 



spin 



mliQl - mi) 



+ 



f-Qif-Qi-^QlQl 



(AA) 



spin 



-I -1 



{M]j-2P-t)^+{mtrt)' 



■{f-b{2f-Pt-P -Mjjf 
{A.h) 



2ReJ2MaMl = -2 



spin 



{Ml - 2P ■ tf + {mtVtf 



Bw{Q\)Bw{Ql) 



(M^ - Ql){Ml -2P-t) + MwTwmtVt 

Qi-f{P-Qif-b-P-bQi-f) 
-Qi-b{P-fQi-f-P-Qif-f) 



2M^f-Pf-b- 



Ql 



(p.fb-f+P-lf-f-P-fb-f) 



+ 



MwTw{Ml -2P-t)+ uitTtiQl - M^) 



{2Qi-f+Qj) 



e'^^^^'PahQiA 



(A.6) 

where g = |e|/sin6'u;, and Bv{Q^) = [{Q^ - M^)"^ + M^T^]~'^ is the vector boson 
Breit-Wigner factor (V = W or Z). For convenience, we take Mz and Tz to be 91 
and 2.5 GeV, respectively, and Mw and T\y to be 80 and 2.1 GeV, respectively. 
The weak mixing angle is taken to be sin^ 9w = 0.23. In all of our rate equations, 
3" is a color factor, with n being the number of quark- antiquark pairs in the final 
state. 
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Note that we have retained the 6-quark mass in the phase space. However, we 
have omitted it in the above matrix elements since to the same order in mi)/Mw 
there is a further graph given by exchanging the b and t quarks in Fig. 1(b), which 
we have not included. We take m^, = 4.7 GeV. For the top width F^, we have 
summed the following expression over (/, /') = {u, d), (c, s), {ue, e"*"), (z^^, {ur, t~^) : 



Fit - - bff') - 7dx. ^^"^^^(3(1 + r)x, - 2x,^ - 4r) 

{t bW btt ) - _ ] _ , ^2^2 + f 2 ^2 



{{l-x, + r)-Ml,Y + Tl,Ml, 

(A.7) 



where xi = 2Ei,/mt,r = mymf, and M^,r^ = M^jz/m^, F^/m^ respectively. 

It is not difficult to analytically integrate out the light fermion pair. The 
resulting expression depends on three variables, which we choose to be the paired 
fermion invariant masses (two of these) and twice the dot product of the Higgs 
four- momentum (P) with the heavy fermion four- momentum (t); we call this 
invariant variable ^ = 2P ■ t . The Lorentz invariant phase space is 



{mt+mhY 



with 



2Qi 



{Mjj+Ql-Ql){Ql+ml-ml)±X'/\Qlmlml)Xy\M]j,QlQl) 



iA.9) 

A technical difficulty arises as this point. As is turned on from zero, ^_ 
decreases, therein going out of the physical region for the = matrix element. 
We have cured this disease by setting ^_ = ^-(m^, = 0),^_|_ = ^-(m^ = 0) -|- 
[^+(m6^0)-e-(m6 7^0)]. 
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We find for the squared matrix element, in the = hmit, 



spin 



+ 
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2Ql{2Ql - - 2(e -Ql- rriDii - - m|) 

{Ml + ml - g2 _ ^)(q2 ^ g2 _ ^ 2g2(g2 _ ^2^ 

(AlO) 



2mf 



m 



spin ^ 



{Ml - + {mtTtf 



{Ql - Qi){2Mfj + mi - 2QI) - M^mf 



+ Ml 



{M-f, + mt){Q-i + mt)-2QlQl 



+ e{Ql-2Qi) 



(A.ll) 



6^2 



m: 



spin 



{Ml - 0' + (mA)' 



(m2, - Ql){Ml - e) + Mw^r^^mA 



• 1$ (^Q? - 2Q2 _ m2 + ^{Ml + QI-Q 



+ 



2M^ 



MlQ\{Ml -Ql- Ql) + m| ((Qf - Qjf - mI{QI + QD) 



+ M^(M^ + m|) + {Ql - Ql){2Ql -Ml- rr?t) 



(A.12) 
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The limit of all massless final state particles is easily obtained from the above 
equations. For the V*V* intermediate state {V = W or Z), the result is 



4 • S^'g^M? 



2Mh 



^ I dlipS^'^'^i^CL (/l • /2) (/{ • +Cr (/i • /^) (/2 • f[) } 

■Bv{2fi-fi)Bv{2f2-f2) 



(A.13) 

For the V = W mode, the chiral couplings are just 

Cl = 1 , Cr = . {A.U) 

For the V = Z mode, they are 

cos^ 9n,CL = ^(^1 + ai)(^2 + ai) + 2^1^20102 , (^.15) 

and 

cos^ 9wCr = ^{vj + al){vl + al) - 2viV2aia2 , (A. 16) 

with Vi = {T^R + TsL - 2Qsm'^0yj)i, ai ^ {T^r - Tsl)^ where T^l and T^r are 
the weak isospin eigenvalues of the left and right helicity fermions, and Q is the 
fermion electric charge in unit of |e|. Integrating these expressions over the two 
polar angles leads to the following expressions for the asymmetry parameters 
(defined in Eqn.(l) in the main text): 





■{M%-Ql-Ql),J^^ , (A17) 
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Mi 



,2\2 



a2iV) 



J dQl\^'\Ml,QlQl)Bv{Q\)Bv{Ql)Q\Ql , 



(yl.l8) 



with 



,2\2 



D = \ j dQl j dQlx'/\Mjj,QlQl)Bv{Ql)Bv{Ql) 



8QlQl+{M]j-Ql-Qlf 
{A.19) 



The final state dependence of the asymmetry ai{Z) is given by 



Cr — Cl 
Cr + Cl 



4viV2aia2 



in + ai){vi + al) 



{A.20) 



In the massless fermion hmit, further integration of the general expression 
(A. 13) may be done analytically, leading to 



Mi 



{Mh-s/Q\ 



,2\2 



dQ\ 



dQl\^/\Ml, Ql Q1)[8Q2q2 ^ ~Ql- Qlf]By{Ql)By{Q 



(A21) 

When kinematics allow one or both vector boson(s) to go on-shell, it will. 
Then the appropriate limits of our general expressions are found by applying the 
NWA, BviQ'^) n5{Q'^ - M^)/MvTv (times a possible symmetry factor of 2, 
as explained in the text). As a check on formulae (A. 17) to (A. 21), we apply 
the NWA to each vector boson and obtain known on-shell results. The on-shell 
values of the asymmetries of ai{V) and a2{V) are given in the text; the on-shell 
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values of Eqn.(A.21) are just 

T{H ^ W+W- ^ /i/{/2/0 = T{H ^ W+W-)BR{W+ ^ fJ[)BR{W- ^ /a/^ 

(A22) 

and 

T{H ^ZZ^ /1/1/2/2) = 2r(ff ^ ZZ)BR{Z ^ fJi)BR{Z ^ /a/a) . 

(yl.23) 

Similarly, the NWA replacement ByiQ^) 7r5{Q'^-M^)/MvTv in Eqns.(A.l- 
6) and (A.8-12) yields the rate for T{H tbW)BR{W ff). 
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FIGURE CAPTIONS 

1. Feynman diagrams for Higgs decay to four fermions. (a) is the gauge 
coupling contribution, (b) is the Yukawa contribution in the case of a massive 
fermion. F is a generic symbol for either W or Z gauge bosons. 

2. Differential decay rate for H — > udsc + udsc versus the scaled W invariant 
mass Minv/Mw, for various scaled Higgs masses Mh/Mw- 

3. Asymmetry parameters ai{V), {i = 1, 2), V — Z, W, for the decays H — > 
Z*Z* — > ddss, and H W*W* udsc versus the Higgs mass Mh- Only ai(Z) 
is dependent on the final state; for uucc {e^e~ ji^ jjT) it is lower than ddss by 
factor of 0.53 (0.029). 

4. T{H ^ tbsc + tbsc) and the massless limit T{H ^ udsc + udsc) versus 
Mjj for constant values of nit/Mu. The latter is roughly 25% of the H WW 
width. The dotted lines correspond to regions excluded by the standard model 
top quark bounds 89 GeV < nit < 250 GeV. The dashed line corresponds to 
T{H bb); we have not incorporated QCD corrections here, so the true bb width 
may be reduced as much as 50% [2] . 

5. T(H tbsc + tbsc) versus mt/Mn, for fixed Mh = 200, 300, 400 GeV. 
The sharp drop in Fjj across the H ti threshold {mt/Mu = 0.5) for increasing 
values of rrit/MH is due to increasing top virtuality. In the range 2mt < Mjj the 
width shown is just equal to 2r{H — > tt)/ BR(t —>■ bsc). 

6. Ti{H — > tbsc + tbsc) versus Mh, where i = Y,G, and Y + G refer to the 
Yukawa, gauge, and total (including the interference term) contributions to the 
partial width respectively; here, mt/MH = 0.6. 
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